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Abstract. We report the XMM-Newton discovery of the first Compton-thick obscured AGN in a Broad Line 
Radio Galaxy, the Gigahertz Peaked-Spectrum source Mkn 668 (OQ-l-208). The remarkably flat 2-10 keV X-ray 
spectrum (observed photon index, F ~ 0.7), alongside with a prominent iron Kq, fluorescent emission line, is a clear 
signature of a Compton-reflection dominated spectrum. Mkn 688 represents a remarkable example of discrepancy 
between X-ray spectral properties and optical classification, as its optical spectrum is characterized by broad and 
asymmetric Balmer lines. The obscuring matter is constrained to be located within the radio hotspots, in turn 
separated by about 10 pc. If the jets are piercing their way through a Compton-thick medium pervading the 
nuclear environment, one could be largely underestimating the radio activity dynamical age determined from the 
observed hotspot recession velocity. The soft X-ray spectrum is dominated by a much steeper component, which 
may be due to nuclear continuum electron scattering, or inverse Compton of the - remarkably large - far infrared 
emission. Soft X-rays are suppressed by a further Compton-thin {Nh ~ 10^^ cm~^) absorbing system, that we 
identify with matter responsible for free-free absorption of the radio lobes. 
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1. Introduction 

Gigahertz Peaked-Spectrum (GPS) radio sources are a 
class of compact (10-100 mas; 10-100 parsecs) radio 
sources, characterized by a simple convex radio spectrum 
peaking around 1 GHz (see O'Dea et al. 1998 for a re- 
view) . Together with the (despite of their name) less com- 
pact and steeper Compact Steep-Spectrum (CSS) sources, 
they constitute a sizable fraction of the 5 GHz selected 
sources, probably as large as 40% IjO'Dea 199811 . 

In the ROSAT AU Sky Survey, GPS/CSS quasars ex- 
hibited a detection rate 3 times lower than radio-loud 
quasars of comparable power IjBaker et al. 1995|l . The ex- 
planation for their X-ray weakness is still matter of de- 
bate. They may represent an intrinsically X-ray weak 
population of radio-loud Active Galactic Nuclei (AGN); 
alternatively, their flux in the soft X-ray band may be 
suppressed by matter heavily obscuring the active nu- 
cleus. The first pointed observations with ROSAT sug- 
gested that X-ray absorption in GPS/CSS quasars may 
be common IjEivis et al. 1 994). The first hard X-ray mea- 



surements of GPS galaxies needed to await ASCA, which 
unveiled column densities of the order of 10^^ cm^^ 
(O'Dea et al. 2000^ .Guainazzi et al. 2 0001. However, the 
number of GPS galaxies for which X-ray spectroscopic 
measurements are available remains to-date dramatically 
low. 

The X-ray properties of GPS galaxies may have im- 
mediate consequences on our understanding of their na- 
ture, and, in turn, on the evolution of large-scale ra- 
dio structures in the universe. The detection of hotspot 
proper motions in an handful of Compact Symmetric 
Objects dPolatidis fc Conway 2003| ), corresponding to dy- 
namical ages ~ lO^-lO'^ years, suggests that at least some 
GPS galaxies should represent an infancy stage of radio 
galaxies proper ([Phillips fc Mutel 1982| ICarvalho 19851 
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IFanti et al. 19951 IReadhead et al. 1996|l . However, there 
is convincing evidence for the presence of large amount 
of matter in the core of GPS sources Ijde Vries et al. 19981 
IPihlstrom et al. 20031 ISnellen et al. 2002|l . It is still un- 
clear whether this matter may have the average density 
(~ 1-10 cm~^) required to ensure the permanent con- 
finement of the radio structure ( |De Young 1993| ), and to 
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prevent ( "frustrate" ) its full development into a Fairnoff- 
Railey (FR) II, or, most likely, into a FR I galaxy 
IjO'Dea et al. 2000|l . X -ray probing of the circumnuclear 
environment may represent an important clue. 

Thanks to the unprecedented collecting area of its op- 
tics, XMM-Newton IjJansen et al. 2001)1 is the ideal ob- 
servatory to perform high-sensitivity, moderate-resolution 
spectroscopy of weak (10~^^ erg cm~^ s~^) X-ray sources. 
An observation program is being carried on, aiming 
at studying the properties of the circumnuclear gas in 
a sizable sample of GPS galaxies. In this paper, we 
report on the observation of Mkn 668, which consti- 
tutes to date one of the very few hard X-ray measure- 
ments of a GPS galaxy ever. XMM-Newton unveiled 
the first, to our knowledge, Gompton-thick AGN in a 
broad-line radio galaxy. Mkn 668 (OQ+208, 1404-f286) 
is one of the closest known GPS galaxy (z = 0.0077; 
[Stanghellini et a l. 1993|). The host galax y shows signs of 
tidal distortion ( |Stanghellini et al. 1993| |, suggesting a re- 
cent merging event. Optically classified as a Seyfert 1 
IjBlake et al. 19701 IMarziani et al. 1998|l . its radio mor- 
phology exhibits two radio lobes along the NE:SW di- 
rection at the approximate distance of 10 pc. Marginal 
evidence for a radio core closer to the former and 
brighter lobe were recently reported by Lister (2003). 
The whole radio structure is embedded in a large diffuse 
halo of diameter ~30 kpc ( |de Bruyn 1990| ). The hotspots 
separate with an apparent velocity of 30-60/ias yr^^ 
UStanghellini et al. 2002 , Li ster 2003, Barone 20 03). This 
identifies Mkn 668 as a potential "young" radio source, 
with an estimated dynamical age of ~100-200 years. 
Kameno et al. (2000) fit the long-wavelength radio cutoff 
with free-free absorption by matter embedding the radio 
structure. They claimed that thermal emission from the 
absorbing gas is responsible for the soft X-ray emission 
detected by ROSAT (Zhang & Marscher 1994). 

An ASGA observation of Mkn 668 
IjGuainazzi et al. 2003|l unveiled a very flat X-ray 
spectrum, together with a possible bright (Equivalent 
Width, EW ~ 900 eV) K^, fluorescent iron fine. These 
results prompted the XMM-Newton observation, which is 
the main topic of this paper. 

In this paper: energies are quoted in the source frame; 
errors are quoted at the 1-ct level for the count rates, 
and at the 90% confldence level for 1 interesting pa- 
rameter for fit parameters and derived quantities; Hq — 
70 km s~i Mpc-^ and go = 0.5. At the Mkn 668 redshift, 
1" ~ 1.5 kpc. 

2. XMM-Newton observation 

XMM-Newton observed the sky region around Mkn 668 
on January 31, 2003. In this paper only results 
from the EPIG ca meras (pn; IStriider et al. 20011 MOS, 
ITurner et al. 200 l|l will be discussed, as Mkn 668 is too 
weak to be detected by the high-resolution spectroscopy 
cameras (RGS). Data were reduced with SAS v5.4.1 
Ij.Tansen et al. 200T)l . using the most updated calibration 



Table 1. Position and count rate (CR) of the sources 
detected by the XMM-Newton EPIC cameras within 3' 
from the optical nucleus of Mrk 668, and having a signal- 
to-noise ratio larger then 5 



0^2000 


^2000 


0.5-10 keV 






CR (lO^^s-^) 


14''07'"00.3"" 


-f28°27'l4"" 


14.5 ±0.4" 


14''06'"57.4'' 


+28°28'5" 


1.41 ±0.19 


14''07'"04.4" 


+28° 29' 37" 


1.03 ±0.18 



"Mkn 668 



files available at the moment the reduction was performed 
(March 2003). Standard procedures were followed for the 
EPIC data reduction and analysis, as detailed in the SAS 
User's Guide HLoiseau 2003|) . In particular: scientific prod- 
ucts were accumulated combining single and double (to 
quadruple) pn (MOS) events (this is most appropriate for 
moderately weak sources to enhance the signal-to-noise 
ratio); background products were extracted from regions 
belonging to the same CCD as Mkn 668 and free form 
contaminating sources; intervals of high particle back- 
ground were removed, applying standard thresholds to the 
E > 10 keV, single event field-of-view light curves (1 and 
0.35 counts per second for the pn and the MOS cameras, 
respectively). After screening, the total exposure time is 
13.2 and 15.5 ks for the pn and MOS cameras, respectively. 
The event lists of the two MOS cameras were merged to- 
gether before extracting any scientific products. Version 
5.1 of the LHEASOFT package was employed for the 
scientific analysis described hereinafter. 

2.1. Imaging 

In the EPIC field-of-view image a source is clearly de- 
tected close to the nominal boresight positions. In the 
pn, 02000 = 14''07'"00.3^ ^2000 = +28°27'14", The 0.5- 
10 keV count rates are (8.6±0.4), and (6.0±0.3)±10-2 s'^ 
in the pn and MOS cameras, respectively. The best-fit cen- 
troid is 1.5" distant from the galaxy core optical position, 
consistent with typical XMM-Newton attitude reconstruc- 
tion accuracy. There is no evidence for intrinsic extension 
beyond the moderately broad XMM-Newton mirror Point 
Spread Function (6" FuU Width Hah Maximum; FWHM). 
Two additional sources are detected at a signal-to-noise 
ratio larger than 5 within 3' from the Mrk 668 nucleus 
{i.e. within the typical ASGA aperture; see Sect. 3.1). 
Their position and count rates are reported in Tab. ^ 
The sum of their count rates amounts at 29 ± 10% and 
< 11% of Mrk 668 count rate in the 0.5-2 and 2-10 keV 
energy bands, respectively. 

2.2. Spectral analysis of the nuclear emission 

Mkn 668 nucleus spectra and light curves were extracted 
from regions of 40" and 25" radius for the pn and the 
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Table 2. Best-fit parameters and results for global fits to the nuclear spectra of Mkn 668. 
Models legenda: #1 = BKNPO+GA; #2 = WA''°/*x(P0'*°/*+2xGA)+WA''"'^'*x(P0^°'^''+GA); 
#3 = WAx(P0"°^*+PEXRAV+3xGA); #4 = WAx (MEKAL+PEXRAV+GA), where: PO = power-law, 
BKNPO — broken power-law; PEXRAV — Compton-reflection; MEKAL = thermal plasma emission; GA = Gaussian 
emission line; WA = photoelectric absorption 



Model 


(102i%i-2) 


A r Aiard 

(10^=* cm- 2) 


Thard 


-fi'brcak 

(keV) 


Psoft 


kT 
(keV) 




XMM-Newton 














#1 


i.3±l:g 




0-7±oi 


2.3±J]:! 


2.i±g:^ 




34.0/41 


#2 


i.2±g:i 


2A±hi 


2.o±g:l 




= Thard 




31.1/41 


#3a* 


1.1 ±0.5 




2.21±8i!^ 




- — Thard 




34.8/42 


#3b 






2V,t 




2.6 ±0.5 




32.2/42 


#4 


o.9±g:i 




2.2 ±0.4^ 






i.3±n 


33.2/42 


ASCA 
















#3a* 


< 0.8 








— Thard 




79.7/77 



^intrinsic, for the Compton-reflection model; * "baseline model" in text; *fixed. 



MOS, respectively. Total source counts are 880 ± 30 and 
620 ± 20 in the pn and MOS, respectively, 15% and 11% 
of which are due to background. No clear evidence for 
variability of the Mkn 668 nucleus flux is seen during the 
XMM-Newton observation. A fit with a constant function 
on the pn 1024 s binned light curves in the 0.5-2 keV and 
2-10 keV energy bands yields xt — ^-^^ ^^"^ xt = 1-01, 
respectively. We will therefore focus in the following on 
the time-averaged spectra only. Spectra were binned in 
order to have at least 25 counts in each spectral channel, 
and to oversample the intrinsic energy resolution of each 
EPIC camera by a factor not larger than 3. Spectral fits 
were performed in the 0.35-15 keV and 0.5-10 keV en- 
ergy ranges for the pn and MOS, respectively A fit with a 
one-component model^ is clearly inadequate (if a power- 
law model is used: — 111-0/47 ly, where v indicates 
the number of degrees of freedom; the residuals against 
this fit are shown in Fig.nj. Phenomenologically, a good 
fit can be obtained with a two continuum model such a 
broken power-law (cf. Tab. provided that a Gaussian 
emission line is included to account for a narrow-band 
feature around 6 keV (observer's frame). The F-test for 
the addition of emission line feature to the broken power- 
law continuum is highly significant: Ax'^/Av — 36.5/3, 
which corresponds to a 99.9992% confidence level accord- 
ing to the F-test^ A zoom of the spectral area around 
6 keV is shown in Fig. [21 (only pn data are shown for 
clarity), where a constant linear rebinning in the energy 
space has been applied to prevent weak, narrow-band fea- 
tures from being visually smeared. The emission feature 
clearly stands out against the underlying continuum. We 

^ All the models discussed in this paper are photoabsorbed 
by cold matter, whose column density is held fixed to the con- 
tribution of our Galaxy along the line-of-sight to Mkn 668: 
iVH,Gai = 1.4 X 10^" cm"^ ^Dickey fc Lockman 1990| l 

^ Following Protassov et al. (2002) we have applied the F- 
test without imposing any constraints on the Gaussian profile 
normalization sign. 
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Fig. 1. Spectra {upper panel) and residuals in units of 
data/model ratio {lower panel) when a photoelectrically 
absorbed power- law model is simultaneously fit to the 
MOS and pn spectra of Mkn 668. Each data point corre- 
sponds to a signal-to-noise ratio larger than 3 

have simultaneously fitted the non background-subtracted 
pn and MOS spectra in the 5.25-7.25 keV energy band, 
using the C-statistics (which is appropriate for spectra, 
whose channel photon distribution is purely Poissonian). 
In this fit, the local pn (MOS) background has been mod- 
eled with a simple power-law (a good approximation in 
such a small energy range), with normalization 1.41 x 10~^ 
(2.61 X lO-'^) keV cm"^ g-^ and photon index 0.38 (- 
0.52). The local continuum has been approximated by a 
power-law as well. A single narrow emission line leaves 
residuals bluewards the energy centroid. They can be ac- 
counted for either by a single broad line, with centroid 
energy = 6.44±[!:g| keV, a = 100±|g° eV and inten- 
sity Jc = (4.2±j!) g) X 10-^ cm-^ s-^, or by a combination 
of three narrow lines, whose energies are well consistent 
with neutral. He-like, and H-like iron. (cf. Tab. OJ. Total 
net counts in the line are 110 ± 50. Although these "local 
fit" results may be affected by uncertainties in the deter- 
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Fig. 2. pn non-background subtracted spectrum in the 
5.25-7.25 keV energy band. Data points have been re- 
binned with a constant width of 75 eV (which corresponds 
to about one-third of the intrinsic energy resolution) . The 
solid line represents a fit with a power-law continuum and 
an unresolved Gaussian profile. The background level is 
indicated by the dashed line 

Table 3. Properties of the emission line components for 
a "local" (power-law and narrow Gaussian profiles in the 
5.25-7.25 keV energy band) and a "global fit" (model #3a 
in Tab. El). Errors are l-cr for one interesting parameter. 
EWs are calculated against the proper continuum (details 
in text) 



"Local fit" 



'Global fit" 



(keV) 


6.42±°;«| 


6.41" 


/<^) (10-« s-i) 


3.4±?:!; 


2.7±S!i 


EW^^'> (eV) 




630 


Ei""^ (keV) 


6.70±E!:}i 


6.7+ 


(10-« cm-^ s-^) 


o.9±o:? 


2.3±U 


EW'-'^'> (eV) 




3400 


Ei^^ (keV) 


6.97±g:gl 


6.96+ 


/P' (10-'' s-i) 


i.2±[!:? 


< 1.0 


EW'-^'' (eV) 




< 1600 



+ fixed 

mination of the local continuum, and need therefore to 
be regarded with caution, they represent a pathfinder to 
guide the global fitting. 

The line complex has a very large Equivalent Width 
(EW = 1.5±g:| keV). Equivalent widths of this or- 
der can be most naturally explained if the continuum 
underlying the emission line is partially or totally ob- 
scured, whereas the line photons reach us through an 
optical path, which does not intercept the obscuring 
matter l|Turner et al. 1997[IMatt et al. 2000|l . In the bro- 
ken power-law fit of the Mkn 668 spectra, the pho- 
ton index above ~2.3 keV (Fhard = 0.7 ± 0.3) is in- 
deed flatter than typically observed in unobscured AGN 
UNandra et al. 19971 IReeves fc Turner 2000|l . and there- 
fore advocates for obscuration of the nuclear emission. 



In light of the above results, we have tested a sce- 
nario whereby the hard X-ray continuum is produced 
by an heavily absorbed nuclear continuum, whereas the 
soft excess and the He- and H-like iron lines are pro- 
duced by electron scattering of the same nuclear contin- 
uum. In practical terms, we have substituted the bro- 
ken power-law continuum with the combination of two 
power-laws, photoabsorbed by column densities N^'^'^ and 
jyhard^ respectively. Even if the soft and hard photon in- 
dices are tied together, the fit is good (x^ — 31.5/42). The 
flat high-energy spectrum is well represented by a "stan- 
dard" AGN spectral index (Ehard = 2.0±g:|) absorbed 
by Nll^'''^ = 2.4 ±J;g xl023 cm-2. In this scenario, how- 
ever, the equivalent width of the neutral iron line {EW ~ 
700 eV) remains unexplained ( Leahy fc Creighton 1993| ). 
A comparatively good fit {x = 35.7/43) is obtained 
if the hard X-ray emission is modeled with a "bare" 
Compton-reflection component (model pexrav in Xspec; 
[Magdziarz fc: Zdziarski 1995| ), therefore assuming that a 
Compton-thick absorber (A^h > cr^^ ~ 1.5 x 10^^ cm"^), 
totally suppresses the nuclear emission below 10 keV 
IIMatt et al. 1999b|l . The large iron line EW 630 eV) is 
naturally produced in this scenario ( Ghi sellini et al. 19941 
IKrolik et al. 1994)l . Any transmitted component is con- 
strained to be absorbed by N^^^''^ ^ 9 x 10^^ cm"^ 
(for r — 2.21 and a warm scattering fraction of 5%). 
In Tab. |31 we summarize the properties of the emission 
line complex in the "global fit" model, if one assumes 
a 3-components (neutral. He-like and H-like iron) de- 
composition. The EW are calculated with respect to the 
continuum against which they are supposedly produced: 
the Compton-reflection continuum for the neutral com- 
ponent, the scattered continuum for the ionized compo- 
nents. The normalizations of the continuum components 
in this model are: iV™fl<=<=ti°n ^ (1.6 ± 0.5) x lO^^^ and 
^scattering ^ (3.5±']:|) X 10"^ kcV cm'^ s'^ , for the 
Compton-reflection and the scattering component, respec- 
tively. They corresponds to observed fluxes of (0.6±q-2) 
and (3.1±J:?) x IQ-^^ erg cm^^ s^^ in the 0.5-2 keV and 
2-10 keV energy bands, respectively, and to unabsorbed 
luminosities of (0.8±g:|) and {■i.l±l:l) x lO^^ cm'^ in the 
same energy bands. We will refer to this model (#3a in 
Tab.^ as "baseline model" hereinafter. 

An alternative origin for the ionized components of 
the emission line complex is thermal emission in an opti- 
cally thin, coUisionally ionized plasma (jBoUer et al. 2 003|l . 
We have therefore tried to fit the spectrum with a com- 
bination of a Compton-reflection continuum, a single nar- 
row Fei emission line, and thermal emission component, 
using for the last the mekal implementation in Xspec 
IjMewe et al. 1985|l . The fit is again statistically accept- 
able. However, the temperature of the thermal component 
(fcT ~ 1.3 keV, with Z < 0.14^©) is not high enough to 
contribute significantly to the ionized lines emission in the 
iron regime. The addition of a further thermal component 
is not required from the statistical point of view. It is 
also worth mentioning that a model constituted by two 
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thermal components does not fit well the EPIC spectra 
{x^/v^ 65.2/42). 

3. Comparison with previous X-ray observations 

3.1. ASCA observations 

Mkn 668 was observed by ASCA between January 4 and 
January 6, 1999. PreUminary results on this observation 
are discussed by Guainazzi et al. (2003), who stress that, 
due to the low statistics, it is impossible to obtain a 
fully unambiguous spectral deconvolution. We extracted 
spectra from the screened event lists available at the 
HEASARC ASCA archive, employing circular regions of 
3' radius around the source centroid, after retaining only 
SIS events with standard (0, 2, 3, 4) grades. Total expo- 
sure time was 65.5 ks (72.1 ks) for the SIS (CIS) cameras. 
Background spectra were extracted from field-of-view re- 
gions, free from contaminating sources. Response matrices 
appropriate for the extracted spectra were generated with 
Lheasoft v5.1. Applying the baseline model simultane- 
ously to the spectra of all the ASCA instruments one gets 
an acceptable fit {x^ /v = 79.7/77). Both the continuum 
and the line best-fit parameters are consistent with those 
derived from the XMM-Newton observation, within the 
admittedly large statistical uncertainties. In particular, a 
large EW iron line is marginally detected in the ASCA 
observation as well, suggesting, alongside with the flat 
hard X-ray continuum, a Compton-reflection dominated 
spectrum. The parameters of the iron line complex, when 
fit with a single Gaussian profile, are: Ec = 6.7±q | keV, 
cr = 250±i[j eV, 7c = (7±t) x IQ-^ cm^^ g-i (correspond- 
ing to an EW = 850 eV). The 2-10 keV flux observed 
by ASCA is almost double as measured during the XMM- 
Newton observation: (6.1±J-2) x 10~^^ erg cm~^ s~^. This 
is suggestive of an historical change of the nuclear emis- 
sion output power by a factor of 2 between the two ob- 
servations, although a contribution from a larger bright- 
ening of one of the sources detected by XMM-Newton 
and encompassed in the large ASCA aperture cannot be 
ruled out. On the other hand, the 0.5-2 keV ASCA flux 
is consistent with the later XMM-Newton measurement: 
(0.9±[;:|) X 10-" erg cm'^ s'^ 

An earlier (July 17 1998) and shorter (CIS exposure 
time 33 ks) ASCA observation produced a 6a detection, 
corresponding to a total 0.5-10 keV unabsorbed flux of 
(5 ±2) X 10"^'^ erg cm~^ s~^. No useful spectral informa- 
tion can be derived from this observation. 



3.2. ROSAT observations 

Mkn 668 was observed 4 times by ROSAT, one of them 
with the Position Sensitive Proportional Counter (PSPC) 
detector, and the remaining with the High Resolution 
Imager (HRI). Preliminary results of the former obser- 
vation are discussed by Zhang & Marscher (1994). We 
have reanalyzed the ROSAT public data, extracting cali- 
brated event lists from the public archive. The results of 



Table 4. Summary of spectral results for the ROSAT ob- 
servations of Mkn 668 (details in text). The superscript 
refers to HRI ("H") or PSPC ("P") observations 



Date 



9/1/92" 
23/1/96" 
10/7/96" 
12/1/1998" 



(10^ 



0.1-2.4 keV 
rate (10"^ s"^) 



2.4:±li 



< 2.4 
2.0 ±0.4 
7.9 ±1.3 
1.6 ±0.3 



1.5 



Ml 
O 



O 



X 



> 



PSPC 



0.5 - 



:992 1995 1998 2001 
Time (years) 



2004 



Fig. 3. Historical 0.5-2 keV flux fight curve of Mkn 668. 
Measurements are corrected for absorption. Error bars 
represent the propagation of the l-cr uncertainties on the 
count rates. Filled circles and upper limit: ROSAT; empty 
squares: ASCA; cross: XMM-Newton 



the analysis are summarized in Tab. 01 The PSPC spec- 
trum can be adequately fit in the 0.1-2.4 keV energy band 
by an absorbed power-law {x^/i' = 17.6/20). In Fig. |21 
we compile eleven years of Mkn 668 soft X-ray flux mea- 
surements. The HRI fluxes have been extrapolated from 
the observed count rate, assuming an absorbed power-law 
with TVh = 1.4 X 10^1 cm-2 and F = 2.2 (best- fit of the 
XMM-Newton EPIC spectra in the 0.5-2 keV energy band 
alone). Caution must be applied in interpreting this plot, 
due to the uncertainties induced by the unknown spectral 
shape in the ROSAT/HRI measurements, and to the po- 
tential contamination of serendipitous sources in the large 
ASCA aperture. Nonetheless, the observed soft X-ray vari- 
ability dynamical range in Mkn 668 is less than a factor 
50%. The XMM-Newton observation has possibly caught 
one of its faintest soft X-ray states ever. 
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Fig. 4. X-ray SEDs for the Circinus Galaxy {z = 0.0015; 
top, pn data) and the Mkn 668 (z = 0.077; bottom, EPIC 
data) nuclei 



of NGC 4051 IGuainazzi et a l. 19981. The X-ray spectrum 
of this Narrow Line Seyfert 1 Galaxy was discovered by 
BeppoSAX to be Compton-reflection dominated during a 
~100 days-long period when the total X-ray output de- 
creased by a factor of ~ 10^ with respect to its standard 
activity level ( |Uttley et al. 1999| ). 

However, no evidence exists in Mkn 668 for large vari- 
ations of the AGN X-ray output. In hard X-rays, only the 
measurements described in this paper are available, which 
limit the variability dynamical range above 2 keV within a 
modest factor of 2. In the soft X-ray range, where a better 
temporal coverage is available, there is no evidence for a 
variability dynamical range larger than 50% (cf . Fig. O • 
There is therefore no direct evidence so far that extreme 
X-ray variability is the culprit for the discrepancy between 
the optical classification and the X-ray spectral properties 
in Mkn 668. 



4. Discussion 

4.1. Mkn 668: a Compton-thick Seyfert 1 galaxy 

The X-ray spectrum of Mkn 668 exhibits two distinct 
components. At energies E ^ 2.5 keV the spectrum is 
remarkably flat. This evidence, alongside with the promi- 
nent Kq, fluorescent emission line, represents a clear in- 
dication for a Compton-reflection dominated spectrum 
IjMatt et al. 20001 and references therein) . The direct view 
of the primary nuclear emission is hindered by photoelec- 
tric absorption with a column density iVn ^ 9 x lO^'^cm"^. 
This is confirmed by the overall similarity between the X- 
ray Spectral Energy Distributions (SEDs) of Mkn 668 and 
of the Circinus Galaxy (cf. Fig.^. The latter is the closest 
Compton-thick AGN known; its flat hard X-ray contin- 
uum and prominent line emission complex are interpreted 
in terms of Compton-reflection of an otherwise invisi- 
ble (below 10 keV) nuclear emission l|Matt et al. 1999;il 
IMoleudi et al.~2008|l. 

If there are little doubts that the X-ray appearance 
of Mkn 668 is that of a highly obscured AGN, its opti- 
cal spectrum shows features more typical of a Broad Line 
Radio Galaxy IjBlake et al. 1970|l . 

Marziani et al. (1993) discuss the peculiar properties 
of the broad optical lines in Mkn 668. The peaks of the 
broad Hq and components are displaced to the red 
with respect to the centroid of the narrow components 
of the same lines. The observed correlation between line 
luminosity and centroid displacement is well explained by 
radiative acceleration in a system of outflowing clouds. 
Marziani et al. (1993) therefore conclude that broad line 
emission in this galaxy partly occurs in a cone of half- 
opening angle ~12° seen at inclination angles i ^ 30°. 

To our knowledge, Mkn 668 is the first known non- 
transient Seyfert 1 galaxy, whose X-ray spectrum is 
Compton-reflection dominated, hence suggesting that the 
nucleus is covered by a Compton-thick absorber. The only 
analogous case reported in the literature is the "off-state" 



4.2. The soft X-ray spectrum 

The soft X-ray spectrum of Mkn 668 is dominated by 
a steeper component, whose origin cannot be unambigu- 
ously unveiled on the statistical basis alone. No clear ev- 
idence for soft X-ray emission lines is found. The 90% 
upper limits on the equivalent widths of mid-Z He-like 
Ka fluorescent Hues {i.e.: Ne ix; Mg xi; Si XIII; S 
xv) are in range 40-80 eV (against the scattering con- 
tinuum), therefore a factor 3 to 10 lower than typi- 
cally observed in "warm scattered" Seyfert 2 galaxies 
l|Guainazzi et al. 1999|l . On the other hand, the profile 
of the emission complex in the iron fluorescent regime, 
significantly more complex than a single unresolved com- 
ponent, suggests a contribution from highly ionized iron 
species. The EW of the iron He-like emission line {EW = 
3.4±^-3 keV) is consistent with optically thin resonant 
scattering (|Matt et al. 1996|l by plasma with standard so- 
lar abundances. A substantial contribution to the ionized 
iron line from thermal plasma is unlikely. 

Although the warm scattering scenario nicely explains 
simultaneously the soft X-ray continuum and the FeXXV 
fluorescent emission, the agreement between the spectral 
energy indices in the radio (aradio — ~1-1) and in the 
soft X-ray (ax — —1-2) may indicate that at least part 
of the soft X-ray photons are produced by non-thermal 
processes, occurring in the same electron population re- 
sponsible for the emission in the radio band. However, if 
we make use of Kameno et al. (2000) best-fit synchrotron 
model for the radio spectrum of the NE lobe, assuming 
free-free absorption with ts = 6.4 ± 0.4, the extrapola- 
tion underpredicts the soft X-ray counts by at least 0.5 
dex (the discrepancy in the hard X-ray band is obviously 
larger). A monitoring of the 4.9 GHz flux in Mkn 668 
between 1975 and 1995 detected only a 20% gradual de- 
crease of the radio flux ( |Stanghelhni et al. 1997| ), ruling 
out an explanation of the discrepancy in terms of lobe 
radio flux variability. There is actually evidence for a 
.steepening of the lobe spectrum with increasing frequency 
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IjDallacasa et al. 2000|) . with aradio decreasing from ~-0.7 
in the 8.4-15 GHz to ~-1.4 in the 15-22 GHz frequency 
range. Finally, the radio spectra available so far do not 
resolve the hotspot, and represent therefore a blending 
between the hotspot and the surrounding lobe. Although 
only strictly simultaneous and aperture-matched observa- 
tions of Mkn 668 in the radio and X-ray bands will be 
able to definitely rule out this possibility, there is no clear 
evidence so far that synchrotron emission by an electron 
population consistent with the 1-15 GHz spectrum repre- 
sents a major contribution to the soft X-ray spectrum. 

Alternatively, Inverse Compton (IC) have been sug- 
gested to be the dominant process in large-scale X- 
ray jets IITavecchio et al. 2nnn i E^lotti et aJ. 2QQ^, par- 
ticularly in GPS quasars ( |Siemiginowska et al. 2002| 
ISiemiginowska et al. 2003| |. The seed photons scat- 
tered off the lobe electrons may be the same syn- 
chrotron photons responsible for the radio hotspot emis- 
sion (Synchrotron Self-Compton, SSC, model). Cosmic 
Microwave Background photons (CMB), or could be lo- 
cally produced in the AGN environment. In the SSC sce- 
nario, the ratio between the IC and the synchrotron lumi- 
nosity is given by the ratio between the photon (itc) and 
the magnetic energy density (um), which in turn can be 
expressed as l|Tavecchio et al. 1998jl : 



6 X IQ-'^i^saT^^ 



(1) 



where Tn is the brightness temperature at the frequency 
1^8.4, in units of 8.4 GHz. At this frequency, Tn = 0.5 for 
Mkn 668 UGhisellini et al. 1993^) . It is therefore implausi- 
ble that this process significantly contributes to the soft 
X-rays. 

The energy density of the CMB at Mkn 668 red- 
shift is ^ 5 X 10~^'^7 erg cm~'^, where 7 is the elec- 
tron Lorcntz factor. One may estimate the magnetic 
field in the hotspot region Bjc applying the Harris & 
Krawczynski (2002) prescription to the observed ratio be- 
tween the radio {S5 ghz — 2.7 Jy; |Stanghelhm et al. 1996') 
and the soft X-ray (Lx = 8 x lO'*^ erg s""'^) luminos- 
ity. Assuming an inclination angle of 15° (45°)"^, one gets 
Bic ^ 3 /iG (5 /iC). This is far too low with respect 
to the minimum energy condition magnetic field (-Bmc; 
[Pacholcyzk 1970] ). The observed 5 GHz integrated fiux 



density in Mkn 668, implies B„ 



60 



-4/7 



(sini)-3A 



( |Miley 1980| ), where 6*860 is the hotspot size in arcseconds. 
In Mkn 668 0sec ~0.7x0.2 mas l|Barone 2003|l . 

The local AGN radiation energy density, uagn poten- 
tially available for Compton scattering should be dom- 
inated by the remarkably large Far Infrared Emission 
(Lfir — 1.4 X 10"L,T^: iMazzarella et al. 199l')l . An exact 
estimate of the energy density at the radio hotspots would 
require a detailed knowledge of the geometry of the cir- 
cumnuclear region, which is largely unknown. If the FIR- 
emitting dust is located within the radio structure (see 

^ Stanghellini et al. (1997) estimate i ~ 45° on the basis of 
the Mkn 668 optical properties, and assuming that the radio 
axis is perpendicular to the plane of the galaxy 



Sect. 4.3.1), and the radio hotspots have a unity filling 
factor we estimate magn ~ 6 x 10^^^(ildust/27r) erg cm~^, 
where fidust is the solid angle subtended by the radio 
hotspots to the dust. For moderate electron Lorentz fac- 
tors, the local AGN radiation energy density can indeed 
exceed the CMB energy density at the hotspots, if dust 
and radio lobes cover regions of comparable size and/or 
are close. 

4.3. The nature of the X-ray absorbers in Mkn 668 

The spectral analysis of the XMM-Newton observation of 
the Mkn 668 unveils the presence of two different absorb- 
ing systems. Together with the Compton thick absorber 
covering the direct line-of-sight to the nucleus, an addi- 
tional absorbing system with Nn ~ 10^^ cm~^ covering 
the soft X-ray spectrum is required by all models pre- 
sented in Sect. 2. In this Section we compare the X-ray 
and radio results, and derive some physical properties of 
the absorbing systems. 

4.3.1. The Compton-thick absorber 

We first ask whether the Compton-thick matter could 
be responsible for the free-free absorption of the ra- 
dio hotspots IjKameno et al. 2000|l . and could there- 
fore be located at a distance from the active nucleus 
at least larger than the half-separation between the 
two lobes (~5 pc). An estimate of the AGN intrin- 
sic power in Mkn 668 can be inferred from the 0[ili] 
fiux (5 X 10"^'' erg cm"^ s~i; IMarziani et al. 1993|) . 
According to the correlation between the 0[iii] and 
the 2-10 keV luminosity [log(L2-io kev/-^o[iii]) ^ 1-8; 
IMulchaey et al. 1994] IMaiolino et al. 1998j . the intrinsic 
2-10 keV AGN luminosity should be ~ 9 x 10^^ erg s"^ 
If the size of the scattering medium is of the same order 
of its distance from the active nucleus, one has: 



4^7 



(2) 



where Isc is the size of the scattering plasma, Lgc is the 
observed 2-10 keV luminosity of the scattering component 
(2 x lO''^ erg s""'^), n is the scatterer particle density, 
^ = Lint/^sc" is the ionization parameter, and il/(47r) 
is the solid angle subtended by the scattering matter to 
the nucleus. If we express: ^ = 500^5oo (the detection of 
He-like iron fluorescent emission line associated with the 
scattering continuum implies a ionization parameter of 
at least a few hundreds); and 17/(47r) = 0.01[ri/(47r)]o.oi 
(corresponding to a cone with a half-opening angle 
of 15°), Isc ~ Coo [f^/(47r)]o.oi 0.02 pc. This would 
be too a small size for the scattered continuum to 
be visible beyond an homogeneous distribution of 
Compton-thick matter covering the radio hotspots and 
the nucleus. This discrepancy is not solved even if 
Q/{4:Tt) = 0.25 is assumed, believed to be a reasonable 
value in the archetypical Compton-thick Seyfert 2 galaxy 
NGC 1068 l|Evans et al. 1991|l . We conclude therefore 
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that the Compton-thick matter must be located within 
the separation between the radio hotspots, in agree- 
ment to typical estimates of the torus size in Seyfert 
galaxies {GrecnhiU et al. 1996", 'Green hill et al. 19971 
2000, Bianchi ct al.. 



2001J . 

From its FIR luminosity, Knapp et al. (1990) estimate 
that ~ 5 X IO^Mq of dust must be present in the core of 
Mkn 668. If this dust is in the form of a torus within 10 
parsecs from the core, the torus column density must ex- 
ceed 3 X 10^^ cm~^ (in agreement with the X-ray results). 
The physical conditions in the torus can be investigated, 
using the stability conditions for X-ray illuminated torii 
discussed by Neufeld et al. (1994) and Maloney (1996). 
The condition for a torus to be fully molecular if that 
log(^off) < —3.4, where <^off is the "effective" ionization 
parameter at the outer side of the torus, taking into ac- 
count self-shielding of the ionizing continuum. Following 
Maloney (1996), ^cff can be expressed as: 



-0.9„-l„-2 



(3) 



44 —1 

erg s , 
Tig is the 



24 



L44 is the X-ray luminosity in units of 10 
is the column density in units of 10^^ ci 
particle density in units of 10^ cm^"^, and rpc is the dis- 
tance between the torus and the active nucleus in parsecs. 
Moreover, a torus is in steady state if its pressure exceeds 
the radiation pressure. If only the gas pressure is relevant: 



ngrp^ > 10 L44Tr^ 



(4) 



where T3 is the average torus temperature (600 K for 
a molecular torus: [Maloney 1996| ). From Eqn01 logrtg + 
21ogrpc > 4.15, which implies log^cff < —4.75, in agree- 
ment with the molecular hypothesis. 

4.3.2. The Compton-thin absorber 

If the jet is opening its way through the ambient in- 
terstellar medium, each radio lobe is surrounded by a 
bow shock, whose precursor highly ionized clouds may 
give rise to free-free absorption of the radio emission. 
Philstrom et al. (2003) estimate from ram pressure ar- 
guments a particle density of ~30 cm"'^ at the exter- 
nal surface of the bow shock generated by the "drilling 
jet" IjBicknell et al. 19 97). If the absorbing matter is dis- 
tributed in a shell surrounding the radio hotspots, its ra- 
dial thickness Arab is 



Arab = — =i 10 pc 



(5) 



This implies that the radio and X-ray (Compton-thin) 
absorbing media could coincide in a layer preceding the 
"drilling" young radio jet. Such a medium would not con- 
tribute significantly to the X-ray emission, even if heated 
up to X-ray emitting temperatures. The expected 0.1- 
10 keV isotropic luminosity of a 1 keV plasma in such 
a geometry would be at most ~ 2 x lO^'' erg s~^. 

The measurement of hotspot recession velocities sug- 
gests that Mkn 668 hosts a jet still in the infancy 




>< 1000 r 



0.1 1 
Compton-thick absorber size (pc) 



Fig. 5. Jet expansion time as a function of the physical 
thickness of the Compton-thick absorber. Solid lines corre- 
spond to I = 15°, dashed- dotted lines to t — 45°. From top 
to bottom, curves plotted in the same style correspond to 
values of the jet opening angle of 20°, 10°, and 5°, respec- 
tively. The light shaded area corresponds to the dynamical 
age estimate from the hotspot expansion velocity. 



of its development ( |Stanghellini et al. 2002|ILister 2003|l . 
Interestingly enough, however, the post-shock particle 
density exceeds the average value (n ~ 1 cm~^), suffi- 
cient to ensure the confinement of the radio jet on scales 
^ 1 kpc ( |De Young 1993| ). The ultimate fate of the jet in 
Mkn 668 may therefore not have been written yet. 

In the light of the above results, we can revise the issue 
of the parsec-scale radio activity age in Mkn 668. Carvalho 
(1998) demonstrated that the presence of a dense medium 
in the core of the host galaxy can significantly delay the 
propagation of the jet, leading to grossly inaccurate esti- 
mates of the evolution time scales. If the jet propagates 
under ram pressure equilibrium through a homogeneous 
external medium, its expansion velocity Vj can de derived 
from HCarvalho 1998|l 



(6) 



where Uo is the electron density and u the jet total internal 
energy density. The expansion time to a distance Rq can 
be therefore expressed as (|Scheuer 19741 ICarvalho 1985|> : 



io - 4 X 10'>R'a ''{cosi)-' ''Lr°-^n°-^NP-^ 



(7) 



where Linj in the luminosity injected in the jet, and Q is 
the solid angle subtended by the jet at Rq. We assume for 
Mkn 668 Rp = 5 pc , and Li„j ~ Lkin 2 x 10^'' erg s'^ 
IjCelotti et al. 1997|l . In Fig. we show tr. as a function 
of the physical size of the Compton-thick absorber, in a 
simple radial geometry where a fraction /thicic^o of the 
optical path along the jet is (actually was) occupied by 
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homogeneously distributed matter with number density 
corresponding to a total column density 5 x 10^*^ cm~^, 
and the remaining fraction (l-/thick)^o by matter corre- 
sponding to a total column density 10^^ cm~^. Large in- 
clination angles, or large thickness of the Compton-thick 
absorber even for moderate inclinations can significantly 
brake the jet, leading to an underestimate of its evolution 
time by one-two orders of magnitudes. 

4.4. Are the AGN in GPS galaxies obscured? 

In the left panel of Fig.Elwe show the 2-10 keV luminosity 
versus the 5 GHz luminosity for the radio-galaxy sample of 
Sambruna et al. (1999) and for the GPS galaxies, for which 
hard X-ray measurements are available so far^. The GPS 
galaxies - which span a range of 5 decades in both radio 
and X-ray luminosity - do not show any deviations from 
the behavior of "standard" radio galaxies. On the other 
hand. Fig. shows the dependence of the X-ray column 
density on redshift^. 3 out of 4 GPS galaxies exhibit an 
absorbing column density larger than 10^^ cm~^, whereas 
only ~30% the Sambruna et al. radio galaxies do. Drawing 
firm conclusions on this issue from the currently available 
data is still impossible. However, they are at least con- 
sistent with GPS galaxies hosting AGN covered on the 
average by larger X-ray obscuration than large-scale ra- 
dio galaxies. The ongoing XMM-Newton program will al- 
low us to enlarge the sample of GPS galaxies with high- 
sensitivity hard X-ray spectroscopic measurements, and 
may therefore provide the final answer to this unsettled 
issue 

5. Conclusions 

Despite of its optical appearance as a Seyfert 1 galaxy, 
the X-ray spectrum of the GPS galaxy Mkn 668 exhibits 
typical features of obscured AGN: 

— above ~2.5 keV, the spectrum is flat and characterized 
by a prominent {EW ~ 600 eV) neutral iron emis- 
sion line. This suggests that Compton-reflection domi- 
nates in this energy range, alongside with fluorescence 
above a heavily absorbed primary nuclear continuum. 

— the column density of the absorber covering the 
primary emission is constrained to be larger than 
9 X lO^'^ cm~^; however values larger than several 
10^^ cm^^ are more likely if the absorber reprocesses 
the AGN emission into the remarkably intense ob- 
served FIR luminosity 

— the bulk of the soft X-rays is due to a steeper (F ~ 2.0— 
2.5) power-law, absorbed by a further column density 
Arg>ft ^ 1021 cm . We associate this component with 
the observed fluorescent Kq, emission from He-like iron. 

NGC 1052, Guainazzi et al. 2000; 13454-125, O'Dea et al. 
2000; Mkn 668, this paper; Q2127-f040. Siemiginowska et al., 
in preparation 

^ Whenever more than one absorbing system is detected, the 
plot in Fig. |S] displays that with the largest column density 



— (warm)-scattering of the nuclear radiation is the most 
straightforward explanation for the soft X-ray emission 
and the Fexxv fluorescent iron line. Nonetheless, 
a contribution to soft X-rays from IR photons up- 
scattered by the hotspot electron population cannot 
be ruled out, although a quantitative estimate would 
require a detailed knowledge of the galaxy core geom- 
etry. 

— we locate the Compton-thick absorber in the inner- 
most region of the galaxy core, within the radio 
hotspots. In order to be visible, the scattering plasma 
needs to be located beyond the visible outer rim of the 
Compton-thick absorbing system (what is commonly 
observed in Seyfert 2 galaxies), probably colocated or 
at least along the same line-of-sight as the unobscured 
BLR (what is instead not usual in Seyfert 2 galaxies, 
where BLRs are hidden from out direct view) 

— we identify the X-ray Compton-thin absorber with 
matter responsible for free-free absorption of the ra- 
dio source. In this scenario, this matter is proba- 
bly distributed in a ~10 pc shell of relatively dense 
{n ~ 30 cm~'^) matter compressed by the jet bow 
shock. 

The discrepancy between the optical classification and 
the X-ray spectral properties in this object does not fit the 
0-th order unification scenario for AGN IjAntonucci 1993L 
IBarthel 1989|l . This might be due to an uncorrelated vari- 
ability between the two bands, driven by large amplitude 
variation in the overall energy output of the AGN, as 
already observed in NGC 4051 IjGuainazzi et al. 1998|l . 
or in the so-called "changing-look" Seyfert 2 galaxies 
l|Matt et al. 2003 bl and references therein). However, there 
is no evidence for large variation of the AGN power in 
Mkn 668. Another possibility is that the BLRs are seen 
in reflected, rather than in direct, light (see, e.g., the dis- 
cussion in Matt et al. 2003a). Given the large EW of the 
broad optical lines IMarziani et al. 1993), this would im- 
ply a large covering factor for the reflector. However, the 
discrepancy may be otherwise related to the peculiar prop- 
erties of this class of radio object. "Young" AGN, triggered 
by recent merging episodes, may have more compact torii, 
which do not effectively cover the direct view of the broad 
line region. Alternatively, we may be observing the AGN 
through the line-of-sight recently opened by the "drilling" 
jet, which frees the view of the broad line regions. In this 
case, a geometry such as a warped disk may still ensure 
simultaneous obscuration of the AGN. It is in principle 
possible as well, that Mkn 668 hosts an "intrinsic" type 
2 AGN, and that the BLRs are produced by the interac- 
tion of gas with the expanding radio lobes. Bicknell et al. 
(1997) make specific predictions on the luminosity of [Oiii] 
and Hq-|-[Nii] optical lines produced in this scenario. The 
observed luminosities in Mkn 668 (~ a few 10^° erg s~^; 
Marziani et al. 1993) are in good agreement with their 
model for reasonable values of the ratio between the lobe 
energy and the radio power (1 Jy at 1.7 GHz; Stanghellini 
et al. 1997). However, the optical lines are probably too 
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Fig. 6. Le]t panel: 5 GHz versus 2-10 keV luminosity for the radio galaxies of the Sanibruna et al. (1999) ASCA sample, 
and the GPS galaxies for which hard X-ray measurements are available. Empty circles: broad-line radio galaxies; empty 
squares: narrow-line radio galaxies; empty triangles: radio galaxies; filled circles: GPS galaxies. Right panel: Nn versus 
redshift distribution for the same sample. 



broad to be produced by such a mechanism. The upper 
limit on the velocity of clouds accelerated by a bow shock 
can be estimated by imposing that the shocked gas cool- 
ing time is lower than the radio source dynamical age 
IjO'Dea et al. 2002|l : 

vci ^ 50i?i/^ip{i^iooV//^ < 60 km s"! 

where Rpc is the distance (in parsecs) from the nucleus to 
the bow shock, nps. loo is the post-shock density in units 
of 100 cm^'^, and Vbs is the bow shock velocity in units 
of the speed of light (~0.1 for Mkn 668; Stanghellini et 
al. 2002). On the other hand, the observed FWHM 
is ~6000 km s~^ ( [Eracleous fc Halpern 1994| ). Moreover, 
a broad (FWHM ~1450 km s~^) Hq component was de- 
tected in polarized flux IjGorbett et al. 1998)l . The inter- 
stellar medium should be distributed according to a rather 
peculiar geometry, for bow shock interactions to produce 
a broad polarized line. 

An interesting suggestion was recently put forward by 
Bellamy et al. (2003), who discovered broad Paa line emis- 
sion and a reddened continuum in the young radio galaxy 
PKS 1549-79. The compactness of the radio emission in 
this object, together with the one-side jet morphology, 
indicate that the radio axis is probably aligned close to 
the line-of-sight. According to the 0-th order AGN uni- 
fication scenarios, the BLR and the AGN in this object 
should be unobscured, contrary to evidence. Bellamy et 
al. propose a scenario, whereby PKS 1549-79 is an AGN 
in the early stages of its evolution. The high obscuration 
should be a transitory phase, which will end as gas and 
dust are dissipated by outflowing gas. The analogies be- 
tween PKS 1549-79 and Mkn 668 are remarkable, both 



exhibiting simultaneously an obscured nucleus, an out- 
flowing BLR, and a small-scale radio jet. However, the 
radio morphology of Mkn 668 exhibits clearly two radio- 
lobes, although their power differs by about one order of 
magnitude UKameno et al. 2000|l . This prevents the line- 
of-sight to the AGN in Mkn 668 from being too close to 
the radio axis. Moreover, the detection of a large-scale 
(~30 kpc; |de Bruyn 1990| ) radio halo may indicate that 
Mkn 668 could be undergoing one of several phases of re- 
current activity. 
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